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ABSTRACT 

Single pulses of pulsar radio emission are modeled as superposition of radiation orig- 
inating from many small subsources that are randomly distributed in the emission 
region. The individual subsources are given an intrinsic finite angular beaming. This 
model can produce fluctuations in intensity and pulse profiles that are similar to the 
microstructure observed in some pulsars. Statistics of the phase resolved flux density 
of a simulated single pulse can be approximated by a lognormal distribution, which is 
in good agreement with observations. 



Key words: 
eral 



Plasmas-polarization-radiation mechanisms: nonthermal-pulsars: gen- 



1 INTRODUCTION 

Observations of single radio pulses from pulsars show a 
high degree of polarization and rapid variation in inten- 
sity, often with spiky structures called the microstructure 
(Craft, Comella & Drake 1968; Hankins 1971). Microstruc- 
tures seem to occur simultaneously over a wide range of 
frequencies (e.g. Rickett, Hankins & Cordes 1975), suggest- 
ing that they are closely associated with intrinsic proper- 
ties at the emission origin. Three types of mechanisms have 
been considered as possible causes for the microstructure: 
(1) temporal modulation of the relevant emission process 
or wave propagation, e.g. fluctuations in radiation as the re- 
sult of modulational instability associated with the coherent 
emission process (Harding & Tademaru 1981; Chian & Ken- 
nel 1983; Weatherall 1998) or as the result of modulation of 
radio waves due to a low frequency wave (Machabeli et al 
2001); (2) beaming effect attributed to individual emission 
beams sweeping across the line of sight (Hankins 1972); (3) 
the combined effect of relativistic beaming and a nonuni- 
form distribution of the emitting plasma along the curved 
field lines (Luo & Melrose 2004). The broadband nature of 
microstructures (Lange et al. 1998; Kramer, Johnston & van 
Straten 2002) appears to favor the second or third mecha- 
nisms, which are intrinsic to the emission origin. The third 
one differs from the first in that the mechanism requires 
inhomogeneities in a longitudinal distribution of the emit- 
ting plasma along the curved field lines. Since nonuniform 
structures produce fluctuations within the individual emis- 
sion beam, the third mechanism can lead to much narrower 
microstructures than predicted by the conventional beaming 
model. 

In both beaming models (mechanisms (2) and (3)) for 
microstructures a nonuniform, nonstationary pair cascade 



above the polar cap (PC) is required. In the widely-dicussed 
PC models the radio emission is generally assumed to be 
produced in the outflowing electron/positron pair plasma 
within the pulsar magnetosphere (e.g. Melrose 2000). A 
rotation-induced electric field above the PC accelerates pri- 
mary electrons (or positrons) to ultrarelativistic energies; 
these particles radiate high energy 7-rays, which decay into 
pairs in the strong pulsar magnetic field, forming an outflow- 
ing pair plasma (Sturrock 1971; Arons & Scharlemann 1979; 
Daugherty & Harding 1982). The assumption of steady pair 
production, which is a major feature of the current mod- 
els, may be oversimplifled since particle acceleration and the 
subsequenct pair cascade are strongly affected by an exter- 
nal current flow (e.g. Lyubarskii 1992; Shibata 1997; Mestel 
1998) and the pair cascade can be nonstationary and oc- 
cur nonuniformly across the PC. Therefore, it is plausible 
that the outflowing pair plasma is highly inhomogeneous in 
space and varies rapidly in time. Owing to the beaming, a 
radiation pattern generated in such inhomogeneous, nonsta- 
tionary plasma that streams relativistically along the curved 
magnetic field lines, must have nonuniform transverse struc- 
tures that sweep across the line of sight. This can give rise 
to both the fiuctuations and the pulse structure. 

In this paper sinlge pulse emission is modeled as ra- 
diation from many subsources randomly distributed in the 
emission region. The fluctuations and pulse structure caused 
by the random distribution are considered. Individual sub- 
sources radiate in the field line direction with a finite an- 
gular spread. The observed single pulse is considered as su- 
perposition of many subsources that naturally provide the 
microstructure. Their random distribution leads to fluctu- 
ations in intensity from pulse to pulse. To model emission 
from multiple sources, one assumes that the radio emisison 
is in the plasma natural modes that can escape the pul- 
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sar magnetosphcrc (e.g. Melrose 2000). The radio emission 
propagates away from the emission region with the polariza- 
tion being that of the natural modes up to the polarization 
limiting region (PLR), beyond which the polarization is no 
longer affected by the plasma and is frozen to its value at 
PLR (Mehose & Stoneham 1977; Barnard & Arons 1986). 
Specifically, single pulses are simulated numerically with the 
polarization properties derived from the local plasma disper- 
sion at the PLR. 

The multiple subsource model is described in details in 
Sec. 2. Numerical simulation of single pulses and the im- 
plications for the interpretation of the microstructure and 
fluctuations in intensity are discussed in Sec. 3. Conclusions 
are given in Sec. 4. 



2 RADIATION FROM RANDOMLY 
DISTRIBUTED SUBSOURCES 

Single pulse emission is modeled as superposition of emission 
from a random distribution of subsources in the omission re- 
gion. A nonstationary pair cascade above the PC produces a 
nonsteady, inhomogeneous pulsar plasma. The correspond- 
ing source is then highly inhomogeneous and nonstationary. 
One may model such source in terms of a distribution of 
multiple subsources. Polarized radio emission can be com- 
pletely described by the Stokes parameters (7, U, Q, V). 
Assuming that these subsources are not phase related, the 
Stokes parameters can be written as a sum of those from 
individual emitters, 



i 
i 

Q = ^^^(7+-7r)cos2xi, 

i 



(1) 



where 7* are the intensities of radiation from the ith source 
in two orthogonal modes ±, Xi is the position angle (PA) 

of a ray originating from the ith subsource in the observer's 
direction, and ^' and ^'^ are the degree of linear and circular 
polarization, given by 



1 - 



2T 



(2) 



1 + T2 ' ^ 1 + T2 ' 
with T the polarization ellipticity of the -|- mode. 

2.1 Relativistic beaming 

In the observer's frame, radiation from subsources is beamed 
in the direction of motion of the source with an angular 
spread AOo = l/Ps, where Fs is the bulk Lorentz fac- 
tor of the source (as shown in Figure 1). Both aberra- 
tion and refraction can change the beaming direction sub- 
stantially (Blaskiewicz, Cordes & Wasserman 1991; Petrova 
2000; Fussell & Luo 2004). The latter effect cam also cause 
the two modes to separate (Melrose & Stoneham 1977). All 




Figure 1. Radio emission from multiple subsources. The radio 
waves are emitted from the subsources at a and 6, which are 
beamed in a narrow cone with a half opening angle of 1 /F^ in the 
field line direction. The polarization is assumed to be decoupled 
from the local plasma and field line curvature at a' and b' . 



these effects are ignored in the following discussion. The an- 
gular distribution of the intensity of the individual subsource 
is written as 



loi exp 



2(1 - cos9kb) 



(3) 



where k is the direction of the line of sight, 7^ is the central 
intensity of the beam and Okt is the propagation angle with 
respect to the field line direction b. 

One may express 9kb in terms of the polar angles of 
k and b given respectively by {9k,(l>k) and (9bi(j>b)- In the 
relativistic limit Fs » 1, it is convenient to use the small 
angle approximation \6k — ^bl ^ 1 and \(f>k — 4>b\ ^- Then, 
the propagation angle in the observer's frame is written as 

Okb « [(6'fc - ebf + sinOk smeb{<l>k - 4>bf] 
2.2 Distribution of subsources 

The distribution of subsources is closely related to that of 
the pair cascade above the PC. The dominant pair creation 
process is single photon decay in the pulsar magnetic field, 
which is the most efficient on the field lines with the small- 
est radius of curvature. Therefore, pair production should 
be peaked on the field lines near the surface boundary sub- 
tended by the last open field lines where the electric field 
is strong and the radius of field line curvature is small (e.g. 
Arons & Scharlemann 1979; Daugherty & Harding 1982). In 
contrast, there are few pairs created near the magnetic pole 
where the field line curvature tends to become infinite large, 
or on the last open field lines on which the parallel electric 
field starts to drop off rapidly to become zero (on the closed 
field lines). 

Similar to the distribution of pair production, one as- 
sumes that subsources have a distribution peaked near the 
last open field lines. Since the pulsar plasma flows relativis- 
tically along open field lines, one may assign each subsource 
a magnetic polar coordinate {r,6,4>), where (9, <A) are the 
polar angles of the field lines on the PC, and r the radial 
distance of the source to the star's center. One assumes that 
the distribution is peaked on the field lines 6c = ecOd, where 
< ec < 1 is the parameter that characterises how close 
the peak is to the outer rim of the PC where the last open 
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field lines intercept the stellar surface. The PC is defined by 
the half-opening angle Od = (n7?o/c)^^^, where n = 2n/P, 
Ro = 10^ cm is the star's radius. The general form of the 
distribution is then given by 



/(e,0)=5(0)exp 



1 fO-Oc\ 
'2 V eOc J 



2p 



(4) 



where < 6 < 6d and < (A < 2it, c = {cj , C2} and p > 1 arc 
three parameters that characterise the general shape of the 
distribution centered at tc^d, with ei characterising the slope 
at ^ < and £2 at ^ > ^c- One has a Gaussian distribution 
in 9 for p = 1. Since pair creation is generally not axially 
symmetric with respect to the magnetic pole, one includes 
the azimuthal dependence g{4>) in (4). 

2.3 Polarization 

The polarization is assumed to follow that of the plasma 
natural modes until the PLR beyond which the wave re- 
tains its polarization (Melrose & Stoneham 1977; Melrose & 
Luo 2004a). Thus, the polarization in (1) is evaluated at the 
PLR. Specifically, the polarization ellipse T and PA Xi are 
written as a function of the distance d < Rlc = c/Q, along 
the propagation path from the source to the PLR, called the 
decoupling distance, where Rlc is the light cylinder radius. 
Similar to the conventional rotating vector model (Radhakr- 
ishnan & Cooke 1969), the PA at the PLR can be determined 
by the orientation of the local field line plane (at the PLR). 
Since wc emphasize the microstructures and fluctuations in 
intensity, to simplify numerical calculation the ellipticity T 
is evaluated in the cold electron-positron plasma with a net 
charge density r] = (n+ — n-)/(n+ + n-) ^ 0, where n± 
are the number density of electrons (positrons). A detailed 
discussion of the ellipticity including the effect of the rel- 
ativistic distribuion is given in Melrose & Luo (2004a, b). 
For u) ^ u!p, where is the plasma frequency, the natural 
modes can be regarded approximately as transverse and T is 
approximately a Lorentz invariant. It is then convenient to 
calculate T in the plasma rest frame (Melrose & Luo 2004a). 

In the low frequency approximation, ui <^ Qe, T changes 
its sign if the propagation angle in the plasma rest frame 
sweeps across 90°, which is Lorentz transformed to 9kb ~ 
l/Fs in the observer's frame (Melrose & Luo 2004a). The 
ellipticity can also change sign if the net charge density 
changes its sign, e.g. it may occur on null surface where 
the magnetic field direction is perpendicular to the rotation 
axis (Figure 2d and 2f). 



3 THE RESULT AND DISCUSSION 

In the simulations one draws a random distribution of sub- 
sources from (4) with p(0) = 1, a radial range of r to r -|- Ar 
and angular ranges of < 6^ < and < (ji < 27r to evaluate 
the Stokes parameters (1). The simulated pulse profiles are 
shown in Figure 2 for a pulsar with the inclination a = n/3, 
the pulse period P = 0.1s. The pulsar magnetic field is as- 
sumed to be a dipole with B = 10^^ G on the PC. For the 
efficiency of the numerical calculation, one considers a small 
radial range Ar <C r, where all radial distances are in units 
of the star's radius Ro- The plasma has a bulk Lorentz fac- 
tor Fs = 100 with a net charge density rj = —0.1. A case 



with sign change, i.e. rj = —0.1 for f2 • b < and r; = 0.1 for 
n • b > 0, where f2 is the spin axis, is also coriisdcrcd. One 
further assumes the emission frequency uj/tOp = 20 in the 
pulsar frame and that all subsources have the same central 
intensity Joi — const. The distribution of subsources is char- 
acterised by p = 1 (Gaussian), tc = 0.8 and ei = £2 = 0.1/£c. 
Because the distribution of subsources along the line of sight 
is random, the total intensity as well as the polarized inten- 
sities fluctuate across the pulse phase. 



3.1 One dominating mode 

As an example one considers the case of one dominating 
mode with a-m = + ) ~ 0.9. All six simulated 

pulses in Figure 2 show depolarization in intensity, though 
the individual subsources are 100% polarized. One may es- 
timate the polarized intensity of the emergent radiation as 
Ip < (2a„i-l)/ ^ 0.8/ with I = The depolarization 

is due to both the effects of mode mixing and cancellation of 
polarization of radiation from sources originating from dif- 
ferent divergent field lines. A spread in the radiation beam 
of the individual subsource also leads to depolarization. The 
latter two effects are especially important if the decoupling 
distance d is short and the line of sight samples a large num- 
ber of subsources originating from the field lines that diverge 
away from the viewing direction (cf. Figure 2c). 

PA variation can be distorted because of the cancella- 
tion effect on the polarization within the emission beam, 
which is strong at the emission origin. An example of scat- 
tering in PA is shown Figure 2c. As one chooses a small d, 
due to the strong cancellation effect on the polarization, one 
has broad scattering in PA. Some pulsars arc seen to have 
a large spead in PA (e.g. Gil & Lyne 1995; Karastergiou et 
al. 2002), which can be explained by the combination of a 
short decoupling distance d and a wider angular beaming of 
the subsource. 

The single pulses in Figure 2 have a relatively large 
CP as the parameters used in the simulation correspond 
to the regime of the aberrated backward circular polarizar 
tion (ABCP) (Melrose & Luo 2004a). In this regime, waves 
propagate backward in the local plasma rest frame, and are 
elliptically polarized. 



3.2 Microstructures 

The simulated pulses show substructures that are similar to 
the microstructrue seen in observations. The profiles shown 
in Figure 2 are the snapshots of the radiation pattern pro- 
duced from a random distribution of subsources. The pulse 
structures are due to the transverse (across the field lines) ef- 
fect only. Because of the relativistic beaming each subsource 
has a natural angular width of l/F,,, which gives rise to a 
typical time scale given by (Cordes 1979) 



27rr, 



(5) 



The recent study of cascade above the PC confirms that 
the secondary pairs have a broad distribution peaked at 



a moderate Lorentz factor around F. 



10 , which is not 



particularly sensitive to the pulse period (Zhang & Hard- 
ing 2000; Hibschman & Arons 2001; Arendt & Eilek 2002). 
Hence, (5) predicts a approximately linear relation with P. 
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Figure 2. Simulated single pulses for a pulsar with P = 0.1, a = tt/3. The three (light, solid and dotted) plots in each case correspond 
respectively to the three intensities I, L = (U'^ + Q^)^/"^ , V . The corresppnding PA variation is shown below each profile, a: am = 0.9, 
d = 200, the viewing angle (relative to the rotation axis) i = 71/3+0.4, and the subsources distributed within the radial range r = 50 — 50.5; 
b: As in a, but with d = 250, i = 7r/3 + 0.3; c: As in a, but with d = 5; d: = 0.9, i = 7r/3 + 0.3, r = 42 - 42.5, d = 250; e: d = 280, 
i = 7r/3 + 0.25, r = 20 — 20.5; f: As in e, but with r = 35 — 35.5. The bulk Lorentz factor of the emitting plasma is Ts = 100. The 
subsources are distributed within the colatitude range < 9^ < 8^ centered at = 0.80^1 with a width A9c = O.lff^- 



Observations of microstructures seem consistent with this 
linear relation (Kramer, Johnston & van Straten 2002). For 
T, = 100, P = 0.1 s, one has = 159 /is. 

The pulse structure can be smeared out substantially 
if the subsources are densely packed in the emission region. 
This may explain why some pulsars do not show any mi- 
crostructure. 

Observations appear to suggest that microstructures 
may have a much narrower width corresponding to a time 
scale much shorter than (5) (e.g. Popov et al 2002). However, 
nanosecond structures have been clearly confirmed only for 
giant pulses (Hankins et al 2003), which are thought to have 
a different origin from normal pulses (Romani & Johnston 
2001). Such rapid variation can be produced by rapid tem- 



poral modulation such as plasma turbulence (Hankins et 
al 2003) or the combined effect of the field line curvature 
and an inhomogeneous distribution of sources along the field 
lines (e.g. Gil 1985; Luo & Melrose 2004). In the latter case, 
due to the relativistic beaming and the field line curvature, 
the observer can only see radiation within a very small tem- 
poral window. For a source moving along the curved field 
lines, with a longitudinal extent ALe < 2Rc/Ts, where Rc 
is the radius of curvature, an observer can only see the ra- 
diation within the time interval given by (Luo & Melrose 
2004) 



1 

''''' ^ "Fs ' 
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Figure 3. Distribution of single pulse flux intensities (the his- 
togram plot). The solid line corresponds to a fit to the lognormal 
distribution. The parameters are as in Figure 2a. 



where ur = c/Rc- If the source is nonuniform along the field 
lines over the length scale less than 2Rc/Ts, the radiation 
must have temporal structures with the time scale given 
by (6). Simulation of such rapid variation in intensity can 
be done in the same way as that leads to Figure 2 except 
that one draws different random distributions of subsources 
within one pulse period. 



3.3 Fluctuations in intensity 

Since the number of subsources beaming into the line of sight 
is random, the relevant intensity varies from pulse to pulse. 
Relative intensities of a simulated pulse can be obtained at 
a particular pulse phase by repeating the simulation a large 
number of times. The distribution of the intensities at the 
pulse phase = 0.1 rad near the pulse peak is shown in Fig- 
ure 3. The distribution can be reasonably well fitted by a 
lognormal distribution, which appears to be in good agree- 
ment with observations of single pulses (e.g. Cairns, John- 
ston & Das 2001; Kramer, Johnston & van Straten 2002). In 
this model, one expects the simulated fluctuating intensity 
to follow roughly a lognormal distribution if the dominant 
contribution to the total intensity is from those subsources 
with beaming at an angle, AOi = (1 — 2cos9kb) <Si Af2o, 
with respect to the line of sight. Since the intensity from the 
contributing subsource is given by (3), the total intensity 
can be written as 

I = ^Jo.expf-^') 

. ..,=p(-i^,..||), (T, 

where Iq = ^ . Joi. Assuming a large number of such sub- 
sources, according to the central limit theory (CLT) one 
expects log/ to approach a lognormal distribution. 

The distribution may deviate from lognormal if the total 
intensity results predominantly from those subsources with 
A0i > Af2o, which may occur at the edge of the pulse 
profile. 



4 CONCLUSIONS 

We consider a model for single pulses, based on the hypoth- 
esis that the radiation is superposition of many subsources 
randomly distributed in the open field line region. The 
subsources radiate within a finite beaming angular range 
~ l/Fs. Numerical simulations are carried out to model the 
pulse profile. The model can predict the basic features of mi- 
crostructures and fluctuations in intensity. Simulated single 
pulses show substructures due to individual emission beams 
from subsources that are randomly distributed along the 
line of sight. In this mulitple subsource model, one simu- 
lates fluctuations in intensity by repeating the simulation 
of a random distribution of subsources. The phase resolved 
intensities near the peak of the simulated single pulse are 
well fitted by a lognormal distribution. Other effects such as 
wave absorption/amplification (e.g. Cairns, Johnston & Das 
2001), scattering in an inhomogeneous intervening plasma 
(Ishimaru 1997) can contribute significantly or even domi- 
nantly to the total fluctuations. However, to include all these 
effects in the model one requires specific models for these 
processes in the pulsar magnetosphere. 

The model predicts a large scattering in PA if the PLR 
is close to the emission origin (i.e. the decouple distance is 
short) or the beaming is wide. The spread in PA is due to 
the spread in angles between the lin eof sight and the range 
of diverging field lines from which radiation can be seen at 
any given time. The model reproduces ABCP proposed by 
Melrose & Luo (2004a). The polarization is modeled in a 
way different from the conventional rotating vector model 
(e.g. Radhakrishnan & Cooke 1969; Blaskiewicz, Cordes & 
Wasserman 1991) in that the polarization is elliptical and 
that it is determined at the PLR, not at the emission point. 
The PLR is characterised by its distance d to the emission 
origin and in this model d is treated as a free parameter but 
constrained by d < Rlc- Further work on models in which 
d can be determined quantitatively is needed. 
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